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We pursue a class of visible axion models where the axion mass is enhanced by strong dynamics 
in a mirrored copy of the Standard Model in the line of the idea put forward by Rubakov. In par¬ 
ticular, we examine the consistency of the models with laboratory, astrophysical, and cosmological 
constraints. As a result, viable parameter regions are found, where the mass of the axion is of 
0(100) MeV or above while the Peccei-Quinn breaking scale is at around 10^~® GeV. 


I. INTRODUCTION 


It is well known that the strong interaction conserves CP-symmetry very well: CP-violating processes in the 
Standard Model so far observed can be explained by the phase of the CKM matrix. The CP-conserving nature of the 
strong interaction is, however, quite puzzling since QCD possesses an intrinsic CP-violating parameter, the 0-angle. 
In fact, the effective 0-angle, OeS = 0 + argdetF„ -I- argdetY^, which sets the magnitude of CP-violation in QCD, 
is constrained to be very small, 0eff 10“^*^, from the null observation of the neutron electric dipole moment 
\dn\ < 2.9 X 10 ^®ecm (90% CL) [3]. Here, denote the up-type and down-type Yukawa matrices, respectively. By 
remembering that the phase of the CKM matrix is of 0(1), the above constraint amounts to unnatural cancellation 
between the intrinsic 0-angle and the 0(1) phase of the Yukawa matrices. 

The most attractive solution to this strong CP-problem is based on the Peccei-Quinn (PQ) symmetry [5]. There, 
the P(l) PQ-symmetry is an almost exact symmetry but broken by the axial anomaly of QCD. After spontaneous 
breaking of the PQ-symmetry, the associated Nambu-Goldstone boson, the axion, obtains a non-vanishing potential 
due to the axial anomaly. Eventually, the effective 0-angle is cancelled by the vacuum expectation value (VEV) of 
the axion at the minimum of the axion potential. 

The original realization of the axion 1313 , however, has been excluded experimentally. There, the axion field is 
embedded in Higgs doublets and the decay constant of the axion, fa, is tied to the electroweak breaking scale, vew, 
i.e. fa — Vew- Then the mass of the axion field is roughly given by. 


m 


a 


fa 


0(100) keV . 


( 1 ) 


Here A and ?7i,r denote the decay constant and the mass of the neutral pion, A — 93 MeV and m,r — 135 MeV. Such 
a light axion with fa — vew has been extensively searched for via the decay of mesons and quarkonia, which ends up 
with a lower limit on the decay constant, fa ^ lOTeV (see e.g. i)- 

Laboratory constraints can be evaded if the PQ-symmetry breaking scale is separated from vew and at a scale much 
higher than vew- For such a large decay constant, however, the axion mass becomes very small and has trouble with 
astrophysics. Eventually, the lower limit of the decay constant is pushed up to fa ^ 10®“^° GeV (see e.g. i)- Based 
on these observations, two classes of models of the invisible axion have been proposed, often called KSVZ [101 E] 
and DSFZ [T^ [T3] axion models, and their phenomenological and astrophysical/cosmological properties have been 
extensively studied (for a review, see e.g. [H])- 

In this paper, we pursue another possibility to evade all the constraints, a heavy axion. For that purpose, we need 
another source of the axion mass than the QCD dynamics, i.e. additional breaking of the PQ-symmetry to the axial 
anomaly of QCD. Such additional breaking, however, cannot be arbitrary since newly added PQ-breaking terms spoil 
the successful cancellation of the effective 0-angle at the minimunm of the axion potential. To resolve the dilemma, 
we follow the idea put forward by Rubakov m, where the QCD dynamics in a copy of the Standard Model (mirrored 
Standard Model sector) pushes up the axion mass^. There, the effective 0-angles in both sectors are aligned with 
each other by a softly broken Z 2 exchange symmetry of the Standard Model and its mirrored copy. Thanks to the 
alignment, the 0-angles in the two sectors are cancelled simultaneously at the minimum of the axion potential. In 


^ See |16II17| for recent works on the heavy axion on the line of the Rubakov idea. However, their models have various unsolved cosmological 
problems. See also discussion in sectionlVl 
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this study, we carefully examine whether the idea can be realized consistently with all the constraints, in particular 
with cosmological ones, by constructing a concrete model. 

The organization of the paper is as follows. In sectionjlTj we introduce a concrete model of the axion with the mir¬ 
rored Standard Model sector. There, we also summarize laboratory and astrophysical constraints on axion parameters. 
In sectioning we discuss cosmological constraints on the axion. In section|IVj we discuss cosmological constraints on 
particles in the mirrored sector. In section|Vj we discuss how to differentiate mass scales in the mirrored sector from 
those in the Standard Model sector without spoiling the PQ solution to the strong CP-problem. The final section is 
devoted to conclusions and discussion. 


II. MIRRORED STANDARD MODEL AND AXION PROPERTIES 

Let us first introduce two copies of the Standard Model each of which has a single Higgs doublet. We name them 
the Standard Model sector and the mirrored sector, respectively. In the following, we put primes on parameters and 
fields in the mirrored sector to distinguish them from those in the Standard Model sector. As mentioned above, we 
assume that dimensionless parameters in both sectors are equal with each other due to a Z 2 symmetry. In particular, 
the effective 0-angles in the two sectors are aligned, 0eff = high energy input scale such as the Planck scale. 

The electroweak scale and the QCD scale in the mirrored sector, on the other hand, can be different from those in 
the Standard Model due to a soft breaking of the symmetry (see section]^. 

To realize the PQ-symmetry, we introduce a gauge singlet complex scalar field (j) which couples to vector-like pairs 
of (anti-)fundamental fermions {ql, Qr) and of SU{3)c and SU{3)c' via 

C = + gH'LQR + h-c. , (2) 

where g denotes a coupling constant. Here, we assume that (f> is even under the Z 2 symmetry. This example is nothing 
but an extension of the KSVZ axion model unmn, and the PQ-charges are assigned to be qLQRi—^), and 

1), respectively.^ After (j) obtains a VEV, {(j)) = fa/V^, the argument of (j) becomes an axion field a with a 
decay constant fa- By integrating out the extra quarks, the axion couples to the Standard model and its mirrored 
copy via 


where and denote the field strengths of the S(7(3)c^ and U{l)P gauge fields.^ As we will discuss at the 
end of this section, we assume that the extra quarks mix with the quarks in the Standard model (and they do so 
similarly in the mirrored sector). Thus, they have non-vanishing U{1)y charges, Qy- Let us remind ourselves that 
the effective 0-angles in the two sectors are aligned with each other even at low energies, so that they are cancelled 
at the minimum of the effective potential of the axion. We will confirm this crucial point in sectionfV] 

After chiral symmetry breaking by QCD and QCD' dynamics, the axion obtains a mass through the mixings to 
the pions in the two sectors. In particular, when the dynamical scale of QCD' is much larger than that of QCD, the 
axion mass is dominated by the contribution of the mirrored sector, i.e. 


rria 


1 + 2:' fa 


(4) 


where z' = /md' denotes the ratio of the up to down quark masses in the mirrored sector. Due to the Z 2 symmetry, 
it should be very close to the one in the Standard Model, i.e. z' z 0 . 56 . In this case, a heavy axion with a 
mass of 0 ( 100 ) MeV can be easily achieved, for example, by taking vew — 10 ^ x vew and Aqcd' ~ 10 ^ x Aqcd for 
/a - lO^GeV. 

In the rest of this section, let us summarize laboratory and astrophysical constraints on the axion parameters, nia 
and fa- The crucial difference of the KSVZ-type axion from the original axion model (and the DFSZ-type axion 


^ As is the case of the KSVZ axion model, so-called the domain wall number is one in our model. Domain walls are unstable and hence 
our model is free from the domain wall problem. In viable parameter regions we discuss in the following, domain walls decay much 
before the axion decouples from the thermal bath. Axions produced by the decay of domain walls are absorbed into the thermal bath 
and do not affect the standard cosmology. 

^ The gauge fields are normalized so that the gauge coupling constants appear in their kinetic terms. 
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FIG. 1: Constraints on the axion parameters. The green (light) shaded region labeled by “SN” denotes the constraint from the 
supernova neutrino bnrst duration. The bine (light) shaded region labeled by “HB” denotes the constraint from the lifetime of 
the horizontal branch stars. The pnrple (dark) shaded region labeled by + a” denotes the constraint from the Kaon 

decay. The red shaded region shows the constraint from the proton beam dump experiment CHARM. Two horizontal lines 
show the constraint from the extra qnark search assuming the Yukawa coupling constant in Eq. (§ to be gr = 1 and g — 0.3, 
respectively. 


model) is that the axion couples to the Standard Model sector and the mirrored sector only through Eq. , and does 
not couple to Standard Model fermions at the tree-level. Due to the lack of direct interactions to Standard Model 
fermions, the main decay mode of the axion is the one into two photons through the effective interaction term, 




327r2 



2(4 + ^) \ 
3(1 + z) J 


^{FF + F'F') , 
Ja 


for rria < 3to^. The decay rate of this mode is given by. 


r 


_ 1 



2{A + z) 
3(1 + z) 


/2 ■ 


(5) 

( 6 ) 


For rua ^ 3m,n., a mode into three pions becomes dominant, and eventually, modes into two gluon jets become 
dominant for a much heavier axion, TOq 3> SruTr- The decay rates should be compared with axion models with direct 
fermion couplings, where decay modes are dominated by the modes into electrons and muons for nia > 2me and 
ruo > 2m^, respectively. 

Accordingly, the laboratory constraints on the axion of this type is quite different from those on models with fermion 
couplings (see e.g. |18j for a compilation of the constraints on the axion-like particle with fermion couplings). For 
0-1 GeV, the most stringent constraint comes from the Br[Ar^ —)• 7r=*= -b nothing] < 7.3 x 10 at 90 % CL [T^. 
By remembering that the decay of the Kaon into the axion is caused by the — a mixing, 

Br[Ar"‘= 71=*= -b a invisible)] ~ e^o_aBr[K=‘" —^7r^+ 7r°] , e^o-a ~ ^ , (7) 

Ja[z + 1 ) 

we obtain a constraint, /q > a few TeV for rUa < 0.1 GeV (see Fig.[^.'* This should be contrasted to axion models with 
fermion couplings where the dominant contribution to the Kaon decay comes from the one-loop Penguin diagrams 


In the figure, we approximate that the size of the E949 detector is about 5 m, and we require the axion to travel longer than 5 m before 
it decays to contribute to Br[it'^ tt^ + invisible], although the lower limit on fa does not depend on the precise size of the detector 
significantly. 
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which leads to a tighter limit on fa, fa ^ 0(10) TeV [TH]. It should be also noted that the axion parameters are not 
constrained by rare decay of quarkonia and _B-mesons into the axion due to the lack of fermion couplings.® 

The axion parameters are also constrained by beam dump experiments. Again, however, the constraints are much 
weaker than the case of axion models with fermion couplings. The most stringent constraint comes from the proton 
beam dump experiment CHARM at CERN [20] . In Fig.[^ we translate the constraint in [5D] onto the KSVZ-type 
axion model (the red shaded region). Here, we exclude the parameters which predict at least three events of the axion 
decay within the decay region ranging in distance from 445 m to 480 m from the beam dump target.® It should be 
noted that the constraints are not applicable for rua ^ 3m,r, since the axion decays immediately after it is produced. 

In the figure, we also show astrophysical constraints on the axion parameters. There, the region labeled by “HB” 
denotes the parameter space in which the lifetime of horizontal branch (HB) stars is shorten by the axion production 
via the Primakoff process [3].^ The region labeled by “SN” denotes the parameter space which reduces the SN 1987A 
neutrino burst duration. In the figure, we follow the discussion in |I8). and in particular, we allow the parameters 
with which the mean free path of the axion is much less than the supernova core size of 10 km. As the figure shows, 
astrophysical constraints allow the axion with a mass above 0.1 MeV for fa — 10^"® GeV. 

For completeness, we also show the constraint on the axion parameters from the search for the extra quarks {q^, qji). 
For fa — 0(1) TeV, the extra quarks obtain their masses from Eq. ([^, and hence, they are within the reach of collider 
experiments. In fact, the production cross section of the extra quarks is much larger than 0(1) fb at the 8 TeV 
Large Hadron Collider (LHC) experiment when they are lighter than 1 TeV. In order for the extra quarks to decay 
immediately, we hereafter assume that extra quarks qL mix with down-type quarks (and they do so similarly in the 
mirrored sector) via 


^ = ^zqLdm + f,[q'Ld'm + h-c. , (8) 

where ^s denote small mass mixing parameters and i is the generation index.® Here, we assume that q^[^ has a 
vanishing PQ-charge, so that the mass mixing is consistent with the PQ-symmetry. Through these mixings, the extra 
quarks mainly decay into H-\-b, Z + b and W + t, where we assume that the mixing with the bottom quark is dominant. 
To date, the 95% CL lower limit on the mass of the extra quarks of this type is 640 GeV set by ATLAS collaboration 
at the 8 TeV running with an integrated luminosity 20.3 fb“^ [22j . In the figure, we show the corresponding exclusion 
limit on fa assuming g = 0.3 and g = 1, respectively. This constraint puts the most stringent limit on fa for a heavy 
axion, rUa > 0(100) MeV.® 


III. COSMOLOGICAL CONSTRAINTS ON THE AXION 

In the previous section, we have discussed laboratory and astrophysical constraints of the KSVZ-type axion model. 
We have found that a rather small decay constant fa = 10®”® GeV is consistent with those constraints for rUa > 
0(0.1) MeV. In this section, we discuss whether such parameter regions are consistent with the Standard Cosmology. 

With a rather small decay constant, fa = 10®”® GeV, the axion is kept in thermal equilibrium in the early universe 
via the effective interactions in Eq. ([^. In particular, the axion does not decouple from the thermal bath of the 
Standard Model sector until the Primakoff process freezes-out. In Fig.[^ we show the freeze-out temperature of the 
Primakoff process Tp given in Emu] by horizontal (blue) dashed lines. The figure shows that Tp is lower than the 
QCD phase transition temperature, Tqcd = 0(100) MeV in most of the parameter region. Therefore, the axion could 
affect the Big-Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB) depending on the mass 
and the lifetime of the axion. In the figure, a darker (blue) solid line corresponds to the parameters which satisfy 
Tp ~ iTia- Above this line, the Primakoff process freezes-out when the axion is still relativistic, i.e. Tp > ma- Below 
this line, on the other hand, the axion is kept in thermal equilibrium even at T < ma- There the axion abundance 
gets suppressed by a Boltzmann factor until the temperature decreases down to Tp. 


® In the decay of quarkonia and S-mesons, the axion appears in the final state through the mixing to and hence, branching ratios 
into the axion are highly suppressed. 

® In our analysis, we assume that the efficiency of the axion signal is independent of the mass of the axion and set it to be 0.5. We check 
that our criterion fairly reproduces the constraint at 90% CL in m when we apply it to axion models with fermion couplings [201. 

^ We extract the excluded region from m- 

^ The newly added mixing mass parameters do not affect the effective angles at the tree level since they do not enter the determinant 
of the mass matrices of quarks. 

^ The extra quarks may have a rather long lifetime as long as they do not cause any cosmological problems, and hence, they can be 
stable inside detectors of collider experiments, such as the LHC. In such cases, the lower limit on fa gets slightly tighter due to the null 
results of stable exotic hadron searches m- 



5 








-5 -4 -3-2-10 1 2 

LogiolmJGeV] 


FIG. 2: Constraints on the axion parameters from cosmological arguments. The dark shaded region denotes the exclusion 
limits in Fig.[^ The horizontal (blue) dashed lines denote the freeze-out temperature of the Primakoff process. The diagonal 
(red) dashed lines denote the recoupling temperature via the two photon interactions. The darker (blue) solid line corresponds 
to the parameters which satisfy Tp — rua- The lighter (red) solid line corresponds to the parameters which satisfy Th^ — rua- 
The shaded region below the line of Tre ~ (and labeled by is excluded by the CMB constraints on A'eff. The 

shaded region above the line of Tpe — lOMeV is shown for awareness of the tension to the BBN (in particular to the D/H 
abundance). 


In the figure, diagonal (red) dashed lines show the recoupling temperature T^e which is defined by 


min 


■ rng 

. ’ T . 


xTgC^SH , 


(9) 


where H denotes the Hubble parameter (see also [M])- The darker (red) solid line corresponds to the parameters 
which satisfy T^e ~ ma-^^ For the parameters above this line, the recoupling temperature is below the axion mass, i.e. 
Tfle < m-a, which means that the axion decays after it becomes non-relativistic. In the parameter region below this 
line, on the other hand, the decay and the inverse decay processes of the axion freeze in at Tre > Wa, which makes 
the photon in the mirrored sector recouple to the thermal bath of the Standard Model. There, the axion density 
decreases exponentially by the Boltzmann factor when the temperature becomes lower than TOq. 

Now, let us discuss constraints on the axion parameters from cosmology. First, let us consider the parameter 
region which satisfies Tre > Wa. In this region, the mirrored photon is kept in the thermal equilibrium with the 
Standard Model sector until the temperature gets lower than the axion mass. Thus, the mirrored photon gives sizable 
contribution to the effective number of relativistic species, N^g unless the axion mass is larger than the QCD phase 
transition. In addition, both the photon and the mirrored photon are slightly warmed up by the in-equilibrium decay 
of the axion. Putting these contributions together, we find that N^g deviates from the Standard Model prediction, 
= 3.046 [25] by 


AN,g CS - iV™ ~ I.l , (Tqcd > ma > T^-dec) , (10) 

AN,g = (11/12)^/^ + ~ 0.83 , {T,-dec >mg> , (II) 

where T^-dec = 0{1) MeV and Te-ann = 0(100) keV denote the neutrino decoupling temperature and the annihilation 
temperature of the electron, respectively. We multiply 8/7 to the mirrored photon contribution to account for the 
difference of the bosonic and fermionic contributions to A^eff- For TOq < T^-dec, we have also taken into account the 


In the figure, we find that Tre > rria is always satisfied for rUa > 3m-^ due to a large decay rate. 
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relative dilution of the neutrino contribution to A^eff due to the axion decayd^ By compared with the constraints on 
A^eff from the CMB observation, A^eff = 3.15 ± 0.23JI26j.^^ we hnd that the parameter space of T^e > rua is excluded 
for TOq <C Tqcd- In Fig.j^ the region labeled by denotes the region excluded by the axion decay contribution 

to the A^eff- 

Next, let us consider the parameter region which satisfies Tfje < ma- As shown in the figure, most of such parameter 
space also satishes Tp > nia, and hence, the axion in this parameter region decouples from the thermal bath when 
the axion is relativistic and decays after it gets non-relativistic, that is, the decay of the axion is of out-of equilibrium. 
In addition, contrary to the case of Tpe > ma-, the mirrored photon does not recouple to the thermal bath of the 
Standard Model sector. With these two differences, the decay of the axion contributes to A^eff differently from the 
previous case. To infer the contribution to A^eff, let us notice the energy density relations at around the decay time 
of the axion, T ~ Tp^,, 


P^+e+{v) — P 7 +e+(j/) (Ffte) + 

P7' — ~^PaiTpie) ■ 


Here, PaiTpe) is given by 


Pa ^ 


C[3] 5*s(FRe) ^ ^3 


g*s{TF) 


'^aTpa , 


( 12 ) 

(13) 

(14) 


where 5*5 denotes the effective massless degrees of freedom contributing to the entropy density in the Standard Model. 
From these relations, we obtain 


AA’eff 

AA^eff 

where k is defined by. 


= AT' 


rSM 

''eff 


8 9*s{TRe) K 

7 2 1 + K 


^ , (Tqcd :S> ma> 7b_dec) 


^SM ^ 8 g,s{TRe) ^ I ^ _ ^SM ^ >ma> Ta-ann) 


(15) 

(16) 


1 30C[3] 1 g*s{TRe) mg 

2 TT^ 5*s(rF) g*{TRe) TRe 

Here, 5 * denotes the effective massless degrees of freedom contributing to the energy density in the Standard Model. 
By numerical calculation, we find that k ;< 0.1 in the parameter space for Tr^ < ma, and hence, the resultant AN^g 
is consistent with the constraint from the CMB observation. 

It should be noted that the axion decay also affects the baryon-to-photon ratio g, which alters the predictions of the 
BBN. In the region of Tpg > ma, the baryon-to-photon ratio measured in the CMB observation, ? 7 cmBj corresponds 
to 77 = (1 -|- before the axion decay. By remembering that the primordial D/H abundance which is highly 

sensitive to 77 is measured precisely, even a slight change of 77 leads to inconsistency between the BBN prediction 
and the measurements of the D/H abundance. To derive precise exclusion limits on the axion parameters, however, 
delicate analysis involving the evolution of the axion energy density along the BBN is required, and it goes beyond 
this paper. Here, instead, we lightly shade the region where the decay of the axion could affect the BBN to note this 
issue. 

So far, we have implicitly assumed that the electron in the mirrored sector decouples from the thermal bath before 
the QCD phase transition of the Standard Model sector. By remembering that there is a Primakoff process between 
the axion and the mirrored electron, the decoupling before the QCD phase transition requires either the axion or the 
mirrored electron must be heavier than Tqcd- In these cases, the annihilation of the mirrored electron does not affect 
the above discussion. 

Let us comment on what happens if both the axion and the mirrored electron masses are below the QCD scale. 
In this case the mirrored electron annihilates into the mirrored photon after the QCD phase transition. Then, for 


The ^A^eff is different from the one given in m because the axion in this model decays into both the photon and the mirrored photon 
evenly, and because the mirrored photon contributes to A^eff- 

If we allow the baryon-to-photon ratio change from the best fit value either, the constraint on A^efF gets slightly weaker Eel, although 
the above deviation, ^A^efF — 1? been still excluded even with such weaker constraints. 
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Tpie < TOa (i.e. out-of equilibrium decay), the resultant mirrored photon contributes to leading to ANgg ~ 2 , 
which contradicts with the CMB observations. For Tfje > rria (i.e. in equilibrium decay), the mirrored photon 
eventually recouples to the thermal bath of the Standard Model sector at the temperature below the QCD scale. In 
this case, the resultant mirrored photon from the mirrored electron annihilation can be redistributed between the two 
sectors. Such a parameter space, however, has been excluded already by the constraints on A^eff as discussed above. 

In summary, we have examined the consistency of the model with cosmology. As a result, we have found that: 

• The axion with Tre > rna and iria ^ 0(100) MeV is excluded by the constraint on A^eff of the mirrored photon 
contribution. 

• The axion with Tfjg < rua and < 0(100) MeV could affect the BBN (the D/H abundance), although delicate 
analysis is required. 

• The axion with uia > 0(100) MeV does not cause cosmological problems. 


IV. COSMOLOGICAL CONSTRAINTS ON THE MIRRORED SECTOR 


In this section, let us discuss cosmological constraints on particles in the mirrored sector. Most of unstable particles 
in the mirrored sector decay very fast. Thus, they cause no cosmological problems. Stable particles, 7 ' e', p' and 
n' could, on the other hand, cause serious cosmological problems unless their abundances are sufficiently suppressed. 
As we have already discussed above, for example, the mirrored electron should decouple from the Standard Model 
sector before the QCD phase transition, since otherwise it increases the mirrored photon abundance. 

First, let us discuss the fate of neutrinos in the mirrored sector. In the Standard Model sector, we assume the 
seesaw mechanism to explain the small neutrino mass m- If the same mechanism works in the mirrored sector, the 
neutrino masses in the mirrored sector, m^/, get enhanced by 


m^i 


'^'ew 


X TOj/ . 


(18) 


As we will discuss in the next section, we mainly consider that jvEW ^ I to make the axion heavy enough, i.e. 
Too ^ 0(100) MeV. Thus, the neutrino masses generated by the seesaw mechanism are much larger than those in 
the Standard Model sector. Eventually, the relic density of the mirrored neutrino exceeds the observed dark matter 
density in most parameter region. In order to evade this problem, the seesaw mechanism should not work in the 
mirrored sector. This can be achieved by turning off spontaneous breaking of the B — L symmetry in the mirrored 
sector (see discussion in the next section) and making the Majorana mass of the right-handed neutrino in the mirrored 
sector vanish. 

Once the seesaw mechanism is turned off in the mirrored sector, neutrinos in the mirrored sector obtain the Dirac 
neutrino mass, 


rrii,' 


/ Mp>m^ 
\ '^‘ew 


1/2 

X VeW ) 


(19) 


which can be much heavier than the pion in the mirrored sector. Here, Mr denotes the mass of the right-handed 
neutrino in the Standard Model sector. With these large masses, mirrored neutrinos immediately decay into a pair 
of the electron and the pion in the mirrored sector, v' ^ e' + tt' . Therefore, neutrinos in the mirrored sector do not 
cause cosmological problems as long as the seesaw mechanism in the mirrored sector is turned off. 

Next, let us consider nucleons in the mirrored sector. Due to their large annihilation cross sections into mirrored 
pions, the abundance of mirrored nucleons is highly suppressed. 





( 20 ) 


The freeze-out temperature of the mirrored photon production via off-shell exchanges of the axion is much higher than the QCD scale 
even for fa — ^ TeV. 

Even if the abundance is lower than the observed dark matter density, there is a hot dark matter constraint, < 10-20 eV EH], 

which amounts to < 10 x v^w- 






One caveat is that the relic mirrored proton becomes dark matter with long-range self-interactions since they couple 
to the massless mirrored photond® The mass density fraction of such dark matter is constrained roughly below 
0(1)% [50]. Thus, as long as the mirrored proton is lighter than 0(1-10) TeV, relic nucleons do not cause cosmological 
problems. It should be noted here that the abundance in Eq. (20) assumes no baryon asymmetry in the mirrored 
sector. This assumption is quite natural if we assume that the leptogenesis |29j explains the baryon asymmetry in the 
Standard Model, since the absence of the seesaw mechanism in the mirrored sector automatically means the absence 
of the baryon asymmetry in the mirrored sector. 

Finally, let us discuss the fate of pions in the mirrored sector. Although pions in the Standard Model sector are 
unstable, the charged pions in the mirrored sector are stable since we have assumed that neutrinos are heavier than 
pions in the mirrored sector. The main annihilation mode of the charged pion is the one into the mirrored photon 
with an annihilation cross section, 


27ra' 


/2 


av ~ 


( 21 ) 


where a' denotes the fine-structure constant in the mirrored sector. Accordingly, the relic abundance is roughly given 
by, 





( 22 ) 


This abundance is sufficiently small as a dark matter component with long-range interactions as long as m,r' 
100 GeV. Therefore, mirrored charged pions in this mass range do not cause cosmological problems. 


V. USE OF SOFTLY BROKEN Za SYMMETRY 

So far, we have treated the QCD scale and the electroweak scale in the mirrored sector as free parameters. In this 
section, we discuss how to achieve those mass scales in the mirrored sector by using a softly broken Z 2 symmetry and 
show that the crucial condition 6eS is maintained even after the breakdown. We also discuss how to differentiate 

the nature of spontaneous breaking oi B — L symmetry in the two sectors. 

Before discussing the origin of the scales, however, let us summarize the relation between these scales and the 
axion mass. In Fig.|^ we show the contour plot of the axion mass for /„ = 10^ GeV. In the figure, we choose 
Aqcd — 400 MeV. In the gray shaded regions, the relative sizes between the dynamical scale and the quark masses 
are different from those in the Standard Model sector. In our analysis, we assume that the chiral phase transition 
happens even when more than three quarks in the mirrored sector are lighter than the dynamical scale. We also 
assume that the extra quarks in Eq. ([^ is heavier than the dynamical scale for simplicity. It should be noted that the 
electroweak symmetry in the mirrored sector is mainly broken by the strong QGD^ dynamics in the parameter space 
with Aqcd' ^ ‘^EW- There, the condensation scale of the mirrored Higgs field also becomes O(Aqcd')) hence, 
the nominal parameter vew is meaningless.^® 

The red shaded (right-upper corner) region is excluded since the mirrored pion mass is larger than 100 GeV. In the 
parameter space with > Aqcd', no pion results from the chiral symmetry breaking in the mirrored sector. There, 
the axion mass is no more given by Eq. (j^ but it scales as Aq^^j^, . It should be also noted that the hadron picture is 
no more reliable in this parameter space, and hence, we need separate discussion on cosmological constraints in the 
heavy quark picture. In the figure, we show the rough exclusion limit where the mass of the quarkonium made of u' 
and d! is larger than I TeV and its relic abundance is expected to exceed 0(1)% of the total dark matter density. 

Here, let us emphasize that we do not need to increase vew' from vew by hand to achieve a viable heavy axion. That 
is, for ma > Tqcd, the mirrored electron is not required to be heavier than Tqcd, and hence, there is no requirement 
to have vew ^ vew- As we have mentioned, however, even if we set the mirrored Higgs mass parameter to be equal 
to that of the Standard Model sector, Aqcd' is eventually required to be larger than vew to obtain ma > Tqcd, 
where the Higgs VEV is of O(Aqcd')- Thus, even if we set vew = '>^ew nominally, the actual electroweak scale 
becomes much larger than vew automatically. 


If U{1)y and U{1)y' gauge bosons have a kinetic mixing with each other, the mirrored proton also has a long-range interaction with 
charged Standard model particles. We assume that the kinetic mixing is negligible. 

In the figure, we do not show the region with Aqqjj/ > fa- Even if we take the potential of </, such that ((f>) < Aqqq/, the condensation 
of the mirrored extra quarks induces (</,} ~ Aqcjj/ through the interaction in Eq. l|^. 

The exclusion regions by either m^i > 100 GeV or m„/ -|- m^i > 1 TeV should be understood as rough estimations. 
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Log^o^^acD 'I^qcd] 


FIG. 3: The contour plot of the axion mass for fa = 10^ GeV. In the gray shaded regions, the relative sizes between the 
dynamical scale and the quark masses are different from those in the Standard Model sector. The red shaded (right-upper 
corner) region is excluded since the mirrored pion mass is larger than 100 GeV. In the blue shaded (upper horizontal) region, 
the masses of the mirror quarks exceed 0(1) TeV. Here, we assume that the seesaw mechanism in the mirrored sector is turned 
off, so that the charged pion (or the corresponding quarkonum) is stable. The (red) dashed line corresponds to the Aqcd' 
which is increased purely by the effect of the larger v^w' ■ 


Now, let us discuss how to differentiate the scales of the two sectors by the soft breaking of Z 2 symmetry. For that 
purpose, let us introduce a spurion field, a{^ 0), which changes its sign under the Z 2 symmetry. Here a has a mass 
dimension one. With the help of the spurion, it is possible to achieve ^ m'jj,{a) and allow them to take almost 

any values. Concretely, we may choose Z 2 invariant parameters, ttiq, mi and c, 

= mg-h miCT-h ccr^ , (23) 

= mg — miCT-h ccr^ , (24) 

so that m|i ^ It should be cautioned here that a cannot be arbitrarily large, since it might appear any 

complex phases of parameters in the two sectors suppressed by the reduced Planck scale, i.e. ct/Mpl with opposite 
signs. In particular, the 0-angles in the two sectors may depend on a by 


£ = 


cr 

MpL 


GG- 


a 

MpL 


G'& 


(25) 


with an 0(1) common coefficient. Therefore, there is an upper limit on the size of the spurion. 


a 

MpL 


< 10 ' 


(26) 


so that too large 0eff does not appear in the Standard Model sector at the minimum of the axion potential. 

Next, let us discuss how to achieve a larger dynamical scale in the mirrored sector. As utilized in [Ml HU to achieve 
a larger axion mass, the dynamical scale of the mirrored sector automatically increases by taking vew ^ vew since 
quarks decouple at higher energy scales than the Standard Model sector. In Fig.j^ we show Aqcd' which is increased 
purely by the larger vew as a (red) dashed line. The figure shows, however, that the axion cannot be heavy enough 
unless Vew ^ 10^ GeV where the quark mass in the mirrored sector exceeds 0(1) TeV. Thus, in order to achieve a 
viable axion mass, rua > 0(100) MeV, we need to increase Aqcd' itself. 
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One may suspect that this kind of “fine-tuning” is problematic. In the low scale theory, it indeed seems unnatural. However, generally 
speaking, mass scales are what should be generated dynamically in a UV theory. Since we do not know the UV theory, we allow tuning 
of mass scales. Otherwise, we must worry about the weak scale itself in the first place, but it is beyond the scope of our paper. 
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With the help of cr, the larger Aqcd' can be easily achieved by introducing extra scalar quarks whose masses again 
depend on cr, i.e. 


iVg 

£ = ^(m?(CT)|%p+m|(CT)|g'p) , (27) 

i=l 

where Nq denotes the number of the extra scalar quarks.^® As a simple example, let us choose Ng so that the beta 
functions of the SU{5) gauge coupling constant vanish when the mass of the scalar quark is negligible. In this case, 
the ratio Aqcd'/^QCD is roughly given by rriq'/rriq^^ 

By the similar token, we can differentiate the nature of spontaneous breaking of the B — L symmetry in the two 
sectors. That is, by again assuming that the mass terms of the B — L breaking fields in the two sectors depend 

on CT, we can easily turn on/off spontaneous symmetry breaking in the two sectors. As the size of a is limited from 
above, so is the B — L breaking scale, < 10^^ GeV in the Standard Model. Fortunately, such B — L breaking 

scale is high enough to allow thermal leptogenesis in the Standard Model sector. 

Before closing this section, let us discuss how largely the 0-angles in the two sectors deviate with each other by the 
soft breaking of the Z 2 symmetry. As we have discussed, the tree-level contributions to the differences of the effective 
0-angle are controlled by the size of a, i.e. 0fg^ = 0{a/Mpi^). Once we take the weak interactions into account, 
however, there are another sources of the CP-violation, the CKM, the MNS, and Majorana phases of the two sectors. 
Since the electroweak scale and the structure of the neutrino masses are differentiated between the two sectors, the 
radiative corrections to the effective 0-angles are different in the two sectors. Fortunately, those differences appear 
at most through 0((a2/7r)^) effects further suppressed by flavor mixings and quark masses [3T1I31], and hence, their 
effects on the 0-angles are highly suppressed. Radiative corrections including the mass mixing parameters in Eq. ([^ 
also contributes to the effective 0-angles. Such contributions are again suppressed by ( 02 / 7 ’')^ and quark mixings, and 
hence, the resultant deviation of the angles are very small. Therefore, the uses of the softly broken Z 2 symmetry do 
not spoil the success of the PQ solution to the strong CP-problem. 


VI. CONCLUSIONS AND DISCUSSION 

In this paper, we have pursued a KSVZ-type axion model where the axion mass is enhanced by a strong dynamics 
in the mirrored Standard Model sector. As we have discussed, the model is consistent with all the constraints when 
the mass of the axion is of 0(100) MeV or above even for a relatively low PQ-breaking scale, 10®"® GeV. We have 
also noticed that turning off the seesaw mechanism in the mirrored sector solves the two problems simultaneously, 
the mirrored neutrino abundance and the too large relic mirrored nucleon mass density. We have also shown that the 
mass scales of two sectors can be differentiated systematically by using a softly broken Z 2 symmetry without spoiling 
the PQ solution to the strong CP-problem. 

One unsatisfactory aspect of this model is that the axion is no more the candidate for dark matter. As an interesting 
alternative, the neutron in the mirrored sector might be a dark matter candidate for 77i„' ~ 100 TeV. Here, we assume 
that rrin' ~ 100 TeV is achieved by a large Aqcd', which induces the electroweak scale in the mirrored sector at around 
the similar scale. Interestingly, in this parameter region, the mass difference between the proton and the neutron in 
the mirrored sector is dominated by QED’ quantum corrections, and hence, the neutron is automatically lighter than 
the proton in the mirrored sector. Therefore, it can be a good dark matter candidate since it does not have long-range 
self-interactions. It should be noted, however, that the charged pion mass in the mirrored sector is expected to be 
0(100) GeV in this parameter range, and hence, this possibility might have a tension with the constraint on the 
mass density fraction of matter with long-range interactions |30| . This tension can be easily solved, for example, by 
assuming that there are only two right-handed neutrinos in each sector, so that one of the left-handed neutrinos in 
each sector become massless [331 Ell ■ With this additional assumption, the charged pion in the mirrored sector decays 
into a charged leptons and a massless neutrino, so that it does not contributes to the dark matter density.^® 


Here, the reason why we introduced “scalar” quark is that they do not contribute to the effective 6-angles, although it may be possible 
to consider extra fermionic colored particles without affecting the effective 6 angles. 

The mass rrigi should be at most of the order of — 10^^ GeV, since the size of cr is constrained to be cr/MpL < 10“^^. 

Otherwise, rriq <C m^/ cannot be achieved by fine-tuning. Accordingly, for Aqqqz/Aqcd — 10^, for example, the mass of extra scalar 
quarks in the Standard Model sector is of 0(10®) GeV or smaller, and hence, they are in the thermal bath after inflation if the reheating 
temperature higher than 10® GeV as required by thermal leptogenesis. In such case, we need separate discussions on how to make these 
particles unstable. 

We will explore more generic possibilities of dark matter candidates in the mirrored sector elsewhere. 
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Throughout this paper, we have assumed that the U{1) PQ-symmetry is an almost exact symmetry of the model 
broken only by the axial anomaly. It is generically believed, however, that global symmetries are not respected at 
least by quantum gravity, and hence, the PQ-symmetry may be explicitly broken by Planck suppressed operators. 




(n-h4)!M] 


PL 




n +4 , 




(n > 0) 


(28) 


with K = 0(1). Such higher dimensional operators leave a non-vanishing effective 0-angle at the minimum of the 
axion potential. 




f 

J ( 


n+2 


2("+2)/2(n-b3)! 


(29) 


For dimension five operators (n = 1), we obtain 


2\0eff 


10"^° X K 


fa 

104 GeV 


/10 GeV 

V Wa 


2 


(30) 


which is consistent with the current upper limit on the effective 0-angle if fa ^ 0(10^ — lO^) for nia = 0(0.1 —10) GeV. 
The stability against possible quantum gravity effects is the merit of axion models with a small decay constant and 
a large axion mass. It is interesting that a non-vanishing effective 0-angle may be observed in near future. 

We note that the small decay constant and the large axion mass is also advantageous when one tries to understand 
the PQ symmetry as an accidental one resulting from other exact gauge symmetries (see [35j and references therein). 

In invisible axion models, where fa > 10® GeV, one must forbid PQ-breaking operators up to dimension-ten in order 
not to induce too large deviation of the effective 0-angle. It is not trivial to obtain such a high quality of the 
accidental PQ symmetry. In our model, as we have discussed above, it is enough to forbid PQ symmetry breaking by 
renormalizable interactions. 

In this paper, we did not copy the PQ symmetry and the PQ breaking field (j). It is also possible that there exist a 
mirrored PQ symmetry and a mirrored PQ breaking, with couplings 


^ = gHLQR + + h.c. . 


(31) 


Assuming that the two PQ symmetries are softly broken down to a single PQ symmetry by the interaction 

C = + h.c. , (32) 

we obtain an axion model with vanishing effective 0-angles.^® If the breaking scale is sufficiently small, there are 
two light axions. 

Let us comment on how visible the present axion model is. Due to a small PQ-breaking scale, the axion may be 
searched for at high intensity low energy collider experiments (see [THI|36] for related axion search). particular, 
the new beam dump experiment at GERN, the SHiP experiment, is expected to cover axion parameter regions with 
a shorter lifetime (and hence a heavier axion) than the GHARM experiments [37] . Another interesting possibility is 
the direct production of the axion and the radial component s at the LHG experiments. In fact, since they couple 
to the gluon rather strongly, they have sizable production cross sections.®® Once they are produced at the LHG, the 
axion immediately decays into jets, while s decays into a pair of axions which subsequently decay into jets. When 
the axion mass is of 0(1) GeV or below, s appears as a two-jet resonance in the 0(1) TeV region, which is difficult 
to be distinguished from QCD background processes (see e.g. [38j). When the axion mass if of 0(10-100) GeV, on 
the other hand, s decays into two axions which can be distinguished from QGD background processes by looking for 
peaks in the dijet invariant mass distributions made by the decay of the axion |39| . As we mentioned in section[TT] the 
production of the extra quarks at the LHG experiments is also an interesting possibility of the present axion model. 


22 For example, the baryon number conservation in the Standard Model is accidentally guaranteed by gauge symmetries of the Standard 
Model. 

2® In order for the PQ solution to work, M'^ must be real. This is guaranteed by the Z 2 symmetry. 

2 ^ Unlike the models discussed in I18II36I . the axion in this model mainly decay into hadrons in most parameter space, and hence, we need 
further study. 

25 For example, for fa — 2 TeV, the production cross sections of the axion and its scalar partner s via the gluon fusion process are 0(100) fb 
and 0(0.1—l)fb for rria — 100 GeV and rus ~ 1 TeV, respectively at the 8—14 TeV LHC. 
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It is also possible to search for particles in the mirrored sector. Here, we just list possible detection methods. At 
least, mirrored particles couple to Standard Model particles through the PQ breaking field (j). The PQ breaking field 
(j) is produced at the LHC experiments and decays not only into Standard Model particles but also into mirrored 
particles, which are invisible for detectors in the LHC experiments. A channel with jet(s) plus missing energy may 
be useful in searching for the invisible decay. The possible kinetic mixing between U(1)y and U{\)y' gauge bosons is 
also interesting.Mirrored particles may be produced in the collision of Standard Model particles through the exchange 
of the gauge bosons. Thermal relics of stable charged particles in the mirrored sector (see Sec. |IV| ) may be detectable 
through CHArged Massive Particle searches. 

Before closing this paper, let us also comment that the above discussion can be easily extended to the model 
consisting of two copies of the Standard Model each of which has two Higgs doublets, so that the PQ-symmetry 
is realized as in the original Peccei-Quinn-Weinberg-Wilczek (PQWW) axion model (see also [H]). There, the two 
sectors share a unique PQ-symmetry through quartic couplings between the two Higgs doublets in the two sectors. 
Contrary to the mirrored KSVZ-type model, the axion decay constant fa is tied to the electroweak scale in the 
mirrored sector, i.e. fa — v'ew- quark mass ratios z^''^ in the two sectors can also differ with each other due to 
the difference of the ratios of the vacuum expectation values of the doublets in each sector, i.e. tan/3 ^ tan/3'. By 
repeating our discussion, we will find that the axion should be again heavier than 0(100) MeV, so that the model is 
safely consistent with laboratory, astrophysical, and cosmological constraints. Thus, we will again need to invoke a 
mechanism which achieves Aqcd' S> Aqcd separately from the size of (see section 0. 

One interesting feature of the mirrored PQWW model is that the 17(1)qed' can be also broken spontaneously 
depending on the mass parameters and quartic couplings of the Higgs doublets in the two sectors. In such case, even 
the “charged” particles in the mirrored sector can be good dark matter candidates. We will explore those possibilities 
elsewhere. 

Another interesting phenomenological difference of the PQWW model is that the main decay mode of the axion is 
not the one into three pions but into a pair of muons even for iria ^ Thus, this type of the axion can be more 

visible at the future beam dump experiments such as the SHiP experiment. Furthermore, it is also possible to detect 
this type of axion by searching for displaced vertices inside the detectors of the LHC experiments made by the axion 
decay [JD] . 
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